One pion production in neutrino reactions: including nonresonant background 
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We investigate neutrino induced one pion production on nucleons. The elementary neutrino- 
nucleon cross section is calculated as the sum of the leading Delta pole diagram and several back- 
ground diagrams obtained within the nonlinear sigma model. This approach does not introduce any 
new adjustable parameters, which allows unambiguous predictions for the observables. Considering 
electroproduction experiments as benchmark, the model is shown to be applicable up to pion-nucleon 
invariant mass W < 1.4 GeV and provides a good accuracy. With respect to the total one pion cross 
section, the model predicts the background at the level of 10% for the ptt"*", 30% for pn'^, and 50% 
for nTr"*" final states. The results are compared with experimental data for various differential cross 
sections. Distributions with respect to muon-nucleon and muon-pion invariant masses are presented 
for the first time. The model describes the data quite well, with the discrepancies being of the same 
order as those between different data sets. 



I. INTRODUCTION 

The interest in one pion production in neutrino- 
nucleus reactions has recently been revived in view of 
the current experimental search for neutrino oscillations. 
The neutrino energy spectra for the ongoing and coming 
long baseline neutrino experiments are typically peaked 
in the GeV region, the region where one pion production 
along with the quasielastic scattering gives a major con- 
tribution. Besides being interesting as a separate chan- 
nel, pion production constitutes a noticeable background 
for various processes: the pion can be absorbed in the nu- 
cleus and thus mimic a quasielastic event, in Cherenkov 
detectors 7r° can mimic the outgoing electron. Thus, a 
precise knowledge of the corresponding cross sections is 
a prerequisite for the proper interpretation of the exper- 
imental data. 

Understanding of one pion production includes two as- 
pects: a proper description of the elementary process on 
nucleon and a proper treatment of the nuclear correction. 
Here we will concentrate on the elementary process. 

In electromagnetic processes, the one pion production 
data, being plotted versus the invariant mass of the out- 
going pion and nucleon, is seen as a series of peaks. This 
picture was a basis for the so-called isobar models, in 
which the intermediate state of the reaction was treated 
as a baryon resonance. The first prominent peak was 
shown to originate mainly from the Delta [i-33(1232)] 
excitation. The second broader peak receives contribu- 
tions from the so-called second resonance region, which 
includes Pii(1440), 1513(1520) and 511(1535) resonances. 
In electroproduction the resonance excitations are known 
to be accompanied by the so-called nonresonant back- 
ground, which can also interfere with the resonance con- 
tribution. Because the theoretical structure of the reso- 



nance contributions is known the modern precise exper- 
iments on meson electroproduction allow the separation 
of these contributions and the extraction of the informa- 
tion related to the resonances only; see, for example, [l| 
for a review. This information can be expressed in the 
form of the quasiexperimental "data points" for the in- 
variant helicity amplitudes that characterize resonances 
and exclude background. 

In neutrino production, the corresponding approach 
is, first, complicated by the fact that the cross section 
contains in addition to the vector current contribution 
also an axial one and a vector-axial interference contri- 
bution so that more resonance properties have to be de- 
termined. Second, any such extraction of such properties 
suffers from the absence of precise, high-statistics data. 
Here the data were obtained mainly in the 1980s in bub- 
ble chamber experiments. The most relevant ones are the 
hydrogen and deuterium data from the Argonne National 
Laboratory(ANL) and the Brookhaven (BNL) National 
Laboratory [2, ISJ which all suffer from low statistics (in 
comparison to the electroproduction data) so that only 
integrated and single-differential cross sections were re- 
ported. An additional experimental problem in these and 
all other neutrino experiments is that one cannot fix the 
neutrino energy, but has to deal with broad band neu- 
trino beams. Thus, we are facing the problem of fixing 
both more complicated background and resonant parts 
from a very restricted set of data. 

Within the phenomenological models, the way out of 
this situation was to presume that in the i^p reaction, i.e. 
in the isospin-3/2 channel, there is no background for the 
A"*""*" production (see, for example, [J-lSl). This was moti- 
vated by, first, the measured tt N invariant mass distribu- 
tion in this channel and, second, by the absence of a nu- 
cleon Born term in this isospin-3/2 channel. Within this 



picture (once the vector form factors of the A production 
are considered to be fixed from electroproduction data), 
one can fit the Delta axial form factors and use them fur- 
ther for other channels. Recent progress in this direction 
was achieved by refitting the vector form factors from the 

l^o-date electroproduction data on helicity amplitudes 
and refitting the axial form factors in the combined 
analysis of the ANL and BNL experiments 10|, [ll | . 

Even if the axial form factors are fitted to describe 
the data for the p7T~^ final state, we have to go beyond 
the isobar concept and include background contributions, 
when considering pir*^ and mr^ final states. The sim- 
plest argument comes from the experimental observation 
that the cross sections for these two final states are ap- 
proximately equal, while the A contribution alone gives 
cr(p7r°)/cr(n7r+) == 2. The calculated cross sections are 
also shown to be lower than the experimental data. In- 
cluding higher resonances, in particular the three isospin- 
1/2 states from the second resonance region, increases the 
cross sections and improves the situation somewhat, but 
does not account for the missing strength. The additional 
contributions required can be introduced within the as- 
sumption a{pTT^)/a{mT^) = 1/2, (the so called "isospin- 
1/2" background) [7|]. A similar philosophy was applied 
recently in [9| , where the vector part of the background 
was extracted from electroproduction, as it is described 
by the MAID group 12| , and then the magnitude of the 
background was fitted to the ANL neutrino data. 

The obvious way beyond this simplest picture is to 
treat the background as a sum of Feynman diagrams with 
a pion and a nucleon in the final state. Progress in this 
direction was achieved by Sato, Uno, and Lee Il3ll . and 
recently by Hernandez, Nieves, and Valverde |14| and 
Barbero, Lopez Castro, and Mariano 15|. 



which include inclusive cross section Q^, and W{Ntt) 



In this article we use the model presented in IJ] and 
apply it to electron and neutrino scattering on nucle- 
ons. We consider various kinematic distributions and an- 
alyze their sensitivity to the background contributions. 
With the model at hand, we are also in the position 
to check the phenomenological treatment, used in the 
Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) trans- 
port model [l6|. This model describes nucleon-, nucleus-, 
pion-, and electron- induced collisions from some hundred 
MeV up to hundreds of GeV within one unified frame- 
work. Recently, neutrino-induced interactions were also 
implemented for the energies up to few GeV with the 
results presented in [QJ, I17h2C1| |. The code is written in 
modular FORTRAN and is available for download as an 



16|. 



open source 

All current neutrino-nucleon investigations concen- 
trate on reproducing a limited number of distributions. 



distribution 0-0, S 0, El, El H, Q. At the same 
time, experimentally available data on distributions on 
the muon-pion and muon- nucleon invariant masses [2|, y] , 
which restrict the dynamics of any model even further, 
are ignored. We concentrate on interactions with nucleon 
targets and aim at reproducing those distributions. In- 
teractions with nuclei will be discussed in a forthcoming 
publication and are not considered in this work. 

The article is organized as follows. First, we give a 
short description of the model used. Then we discuss the 
electron interactions in the Delta resonance region, pro- 
ducing p7r° and mr^ final states. We especially consider 
the resonance-background interference and possibility to 
introduce a simplified description of the background simi- 
lar to the one discussed earlier [9| . Afterward, results for 
neutrino interactions are presented and compared with 
the available experimental data. 

II. SHORT DESCRIPTION OF THE MODEL 

In this section we give details about the model we use 
for the lepton-nucleon interactions. We are studying the 
process of one pion production in lepton interactions with 
nucleons, i.e. 

l{k,,)N{p^) ^ l{k\,)Nip'^)TTip^) 

for various isospin final states. 

The 5-fold cross section for one pion production is 
given by 



da 
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dE' d cos OdETT d cos O-ndcf)-^ 4 yJ(pkY 
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The dynamics of the interaction is encoded in the matrix 
elements for electromagnetic (EM) and charged current 
(CC) interactions 



\M^\ = CealccL^'H^,, H^, = j^jl , 



where 



Cem — 



AiraQED 
Q2 ' 



Ccc = -;- cos 



The calculation of the leptonic tensor L^'^ is straight- 
forward and gives the standard result. The hadronic ten- 
sor Hfj^i, refiects the essence of the process taken into ac- 
count. The hadronic current j,y varies from model to 
model and can include various contributions. 



Several authors have proposed to describe the current 
as a coherent sum of several diagrams 13|, |lj, |23| : Delta 
pole (Dp), crossed Delta pole (cDp), nucleon pole (Np), 
crossed nucleon pole (cNp), contact term (CT), pion pole 
(pp), and pion in flight (pF), 



"" (p + q)2 - Mjj + i£ 



.75 X 



[Vl^-A%]u{p) 



(4) 



j = JDp + jcDp + JNp + jcNp + JCT + jpp + jpF 

The diagrams considered are shown in Fig. [1] 

The progress in understanding the background can, 

however, be only achieved, if the new vertices introduced 

are considered as known and do not include adjustable 

parameters. 

Hernandez et al. 



14| (from now on called the HNV 



model) have proposed to use the vertices predicted by an 
effective Lagrangian of the SU{2) nonlinear cr-model. 
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(8) 



Here gA = 1-26 is the axial nucleon coupling and 
/tt = 0.093 GeV is the pion weak decay constant, which 
enter the Lagrangian of the a model. The currents de- 
fined in (I2|) — ^ can be used for electromagnetic and 
weak processes, provided that the corresponding form 
factors and isospin coefficients C^V:<^Ov-Np,cNp.,CT,pp,pF 

are given. In the following we summarize them. 
1. Delta resonance 



FIG. 1: Diagrams representing the A pole and background 
contributions to the one pion production in weak charged cur- 
rent scattering on the nucleon |14|] . 

The details of the model and the amplitudes of the 
diagrams are given in [14], and we repeat them here only 
for convenience: 



J Dp 
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xu{p)Sap{p + q)Vl^^^{p, q)u{p) 



In the HNV model the vertices of A production and 
decay (which enter the Delta pole and crossed Delta pole 
diagrams) are treated on the same theoretical grounds as 
earlier in @, [ill, [22 1 . 

The A production vertex, F'^'^, is described as a vertex 
for the isospin-3/2 resonance production 



3/2+ 



V: 
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75 • 



(9) 



(2) and used in Eqs. ([2]), ([3]). In terms of the phenomenolog- 
ical form factors, the vector part is given by 
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N 



and the axial part by 
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3/2 






q^r) 



ct 



^{9"'q-P 
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(11) 



The caUigraphic CY stands either for the electromagnetic 
transition form factors C/^ with N — p,n or the CC vec- 
tor form factors CY . For the A resonance, they coincide 
and — in line with [14| — we use the fit of [7| : 



cYiQ') 



c(^-"\q2) 



1 



QV4M2' 






1 



2.13 

-1.51 
0.48 



DviQ"^) 1 + QV0.776M 2 ■ 
(1 + Q'^/M'^Y denotes the 



1 + 1 



V4M2' 
1 



The function Dv{Q^) = \i- ^ ^ /iviyj 
dipole function with the vector mass parameter My = 
0.84 GeV. (The axial form factors are relevant only for 
CC interactions.) Notice, that in general, the currents 
for different isospin channels differ from one another by 
Clebsch-Gordon coefficients. In the present work these 
are included in Eqs. ^ — ([5]) and defined in Table H] 
Thus, the form factors are the same for different final 
states. 

The axial form factors are taken to be the same as in 
14| . where C^{Q'^) was fitted to the ANL cross section: 

0.867 1 



C-(Q^) 



Da{Q^ 
1 



1 + Q73M1 



ctiQ'l 



c^m 



0, 



CtiQ^) 



0, 



with DAiQ"^) = (1 + Q'^/Mlf and Ma = 0.985 GeV. 

The value C^(0) = 0.867 obtained in Q is in 
contradiction with the predictions of the off-diagonal 
Goldberger-Treiman relation, which expresses C^ (0) via 
the AA^TT coupling constant /* and gives the value 1.2. 
This relation is based on the partial conservation of axial 
current (PCAC) hypothesis, which was tested in several 
experiments and is shown to be satisfied with an accu- 
racy not worse than 10% 
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26 



When considering both 
ANL and BNL data, the recent fit [lH gives C^(0) = 1, 
which is closer to the PCAC prediction, but still is out- 
side the 10% deviation from it. In [9| the ANL data were 
described with a fit for C^, that fulfiUs PCAC. 

The spin-3/2 projector is taken in conventional Rarita- 
Schwinger form 
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TABLE I: Clebsch-Gordon coefficients for various final states 
in EM and CC interactions 
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and used in Eqs. ©, ©. 

The coupling /* of the AA^tt vertex in jup and jcDp 
currents is determined from the free decay width of the 
A resonance, /* 



1.15 



14 



2. Nucleon 

The vertices with nucleons and pious are described 
within the SU{2) nonlinear cr-model. Within this model 
all the vertices are pointlike, including the coupling 
VNN of a nucleon to the vector or axial current. 



The authors of |14l | choose to introduce nucleon form 
factors to the VNN vertex in a phenomenological way. 
All these form factors are considered to be known and can 
be taken from one of the conventional parameterizations. 

For the nucleon vertex we adopt the standard hadronic 



current V'^ 



4^ 



v^ 

^N 



FiY + 1 



2Mn 



cr"^ qa, 



and the axial part 



■A 



Fp 



_ FAl^-f5 + -T^q^l5 

Mn 



(12) 



(13) 



Here, TY (* = ^t"^) stands either for the EM nucleon 
form factors F/^ with N = p,n or the CC form fac- 
tors FY — Ff — Fp. The electromagnetic Dirac and 
Pauli form factors Ff' can be rewritten in terms of Sachs 
form factors, for which we take the updated BBBA-2007 
parametrization |27| . The axial form factors are relevant 



only for CC reactions and for them we assume a standard 
dipole form with the axial mass determined in 28|: 



2mjf 



Q' 



V 

Fa{Q^) 



FA{Q^)^gA[l + ^] , AfA = 0.999 GeV 
Fp{Q^) 



(14) 



3. Other diagrams 

As soon as nuclcon form factors are introduced in the 
model, the conservation of the vector current 

q'^iJDp+jcDp+JNp + jcNp + jCT + jpp+jpp),^ = (15) 



is no longer fulfilled. The way to compensate in Eq. p5|) 
the nonvanishing terms stemming from JNp and jcNp is 
to introduce the corresponding form factors to jcT and 
jpF- As outlined in [1J|, for charged current processes, 
the corresponding form factors are the same as the weak 
vector nucleon ones. It can be shown, that for electro- 
production processes, they also stay the same, i.e.. 



F5'^(Q2) ^ FS^iQ^) - Ff (02) _ J.n(Q2) ^ 



F;?{Q^) = Jf/(0') = FfiQ^) - F^iQ 



(16) 



Another phenomenological factor, introduced in the 
HNV model, is Fp, which accounts for the p- meson dom- 
inance in the axial parts of the CT and pp currents. 

For CC reactions Fp is given by 



F^^'iQ') = 



1 



1 + Qyml ' 
and for EM interactions it is zero. 



rrip = 0.77 GeV 



F^iQ' 







To summarize, the HNV model 1^ phenomenologi- 
cally extends the "pure" nonlinear SU{2) model, but 
introduces no adjustable parameters. Besides the nu- 
cleon and the pion, the model contains only the Delta 
resonance and thus is applicable to the region below and 
slightly above the Delta peak. 



III. ELECTROPRODUCTION AS 
BENCHMARK FOR NEUTRINOPRODUCTION 

In this section, we present the differential cross section 
results for electron scattering, with the purpose to check 
the accuracy of the model and the range of its applica- 
bility. 



We consider electrons of energy E^ — 1.884 GeV scat- 
tered on protons over the angle 9 — 47.94° {cos 9 — 
0.67) and calculate the double differential cross section 
da/driedEe. 
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FIG. 2: Differential cross section for ep — >■ eA at E^ — 
1.884 GeV, 9^ = 47.94° (a) and at E^ = 2.238 GeV, Oe = 
31.93° compared to the inclusive ep — >■ eX data [29|. 

Fig. [2^ shows the full model cross section (solid line), 
obtained as a coherent sum of all the diagrams. It is 
compared with JLab electroproduction data [29|. No- 
tice, that these data are for the inclusive cross section, 
while our curve is for the one pion production, that is for 
the sum of the pn^ and nn^ final states, only. Below and 
at the A peak, our calculations agree perfectly with the 
data. At the same time, the Delta pole diagram alone 
(dash-dotted line) is noticeably below the data. Above 
the A peak, as expected, the data lie above our curve, 
because other resonances contribute in this region as well 
as other channels (for example, two pion production and 
eta production) become kinematically allowed. The inco- 
herent sum of all the diagrams (short-dashed line) is also 
shown in Fig. [5] Below the A peak the interference effects 
are small, while above the Delta peak the interference is 



strong and negative. 

The similar picture is also shown in Fig. [2]3 for Ee — 
2.238 GeV and 9 = 31.93° (0086* = 0.8487). The agree- 
ment with the data is again very good below and at the 
A peak. 

Conventionally, experimental results for one pion pro- 
duction are shown in the form of the cross section for 
virtual photons 

1 da 



T, dfl'dE' 



(7T + £(JL , 



(17) 



where Ff is the flux of the virtual photon field 



Tt = 



OIQED E' W^ 



'JV 
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27r2 Ee 2mArQ2 i_£' 



and e is the degree of transverse polarization of the pho- 
ton 



1 + 2 1 



Q' 



tan 



Here 6 is the electron scattering angle, v the energy 
transfer, Q^ the squared momentum transfer, and W^ — 
rnj^ + 2mNV — Q^ the invariant mass of the final nucleon- 
pion state. 

Our results are compared to the DESY electron-proton 
scattering data ^] in Fig.|3]for p7r° and mr'^ final states 
(middle and lower panels). Two data sets are available: 

(1) for electron energy E^ — 2.7 GeV and the scattering 
angle of 14°, at the A peak position Q^ = 0.35 GeV^; and 

(2) for electron energy E^ — 3.2 GeV and the scattering 
angle of 21°, at the A peak position Q^ = 1.0 GeV^. 
The figure shows the cross section ([TT]) versus the outgo- 
ing electron energy E' in the A region. Here the higher 
E' values corresponds to lower invariant masses W. We 
present the full model calculations (solid lines) , as well as 
the contribution of the Delta pole alone (dashed lines). 

The predictions of the MAID model 12|, |31| are also 
shown as crosses, which provides us an overall compar- 
ison with modern electroproduction data. The MAID 
model, developed by Mainz theory group, is a state-of- 
the-art unitary isobar model for pion photo- and electro- 
production on the nucleon, which fits more than 70000 
data points on 5 and more fold differential cross sections. 
For the GiBUU code, it provides the resonance ampli- 
tudes for electroproduction. The results of the MAID 
model for double differential cross section can be consid- 
ered as being equivalent to the data. They can, therefore, 
serve as a benchmark for our calculation that — con- 
trary to MAID — contains a theoretically well founded 
description of the background amplitudes. 

The full model calculations show an excellent agree- 
ment with the data [30| for the n7r+ final state. At 



high £", corresponding to the invariant mass region be- 
low the Delta resonance, we observe a noticeable increase 
of the cross section in comparison with the Delta pole di- 
agram, which significantly improves the agreement with 
the data. With decreasing £", the invariant mass W in- 
creases, reaching W — 1.29 GeV at the left end of the 
data points for E^ = 2.7 GeV and W = 1.35 GeV for 
Ee — 3.2 GeV. For the pvr" final state, our full model, as 
well as the MAID model, shows a reasonable agreement 
with the data. 

The data are also available for the sum of the final 
states pyr" + nir^ (triangles in the upper panel of Fig. [3]) , 
they agree with our curves up to W — 1.4 GeV. In 
all cases the full model calculations are very close to the 
MAID results. Thus, the model provides the same level of 
accuracy as the MAID model, which ensures the applica- 
bility of the HNV model to the leptoproduction processes 
at least up to VF < 1.4 GeV. 



IV. NEUTRINOPRODUCTION 

In this section we present our results for neutrinos. 

The double differential cross sections da /dE^d cos 9^ 
for the charged current neutrino reactions versus the 
nucleon-pion invariant mass W are presented in Fig. |4] 
for the incoming neutrino energy Ey — 1 GeV and the 
muon scattering angle cos 9^ = 0.6. As we already men- 
tioned, the form factors used are taken to be same as 
in [14], in particularly we use the same C^. Fig. \^ 
shows the contribution of each diagram to the cross sec- 
tion for the p7r+ final state. For this channel, the cross 
sections for the background diagrams are indeed small in 
comparison with the Delta pole contribution. Fig. |41d 
compares the full model calculation (solid line) with the 
Delta pole only (dash-dotted line) and with the incoher- 
ent sum (short dashed line) of all diagrams. For the 
kinematics considered, the interference effect is negative 
above the A peak, and positive below the A peak. 

Figs. |3J;-f show the same cross section for a neutron 
target for the two possible final states, pir^ and mr^. 
As one can see, the background terms are noticeable in 
these cases. The most important contribution is given 
by the CT diagram, which provides a rather steep rise 
of the cross section at low W. The Np diagram domi- 
nates the background for the nn^ channel and, together 
with the cDp, gives a large contribution for the pTr" one. 
The way the background appears in reactions with a neu- 
tron target can partly be traced to pure isospin relations. 
For example, as can easily be deduced from Table U the 
leading isobar contribution. Dp, for the nTr"*" channel is 
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FIG. 3: Cross section (|17|l for ep — >■ epTr" (middle panel) and ep — >■ eriTr"*" (lower panel) at i?e = 2.7 GeV, ^e = 14° (left 
panel) and E^ = 3.2 GeV, de = 21° (right panel) as a function of outgoing electron energy E' compared to the data [30|. The 
predictions of the MAID model are shown as crosses. 



9 times smaller than that for the p7r+ one. At the same 
time, the cDp term, which was very small for the p7r+, is 
9 times bigger for the mr'^ and thus becomes noticeable. 
^From Figs. |Hi,f one can see, that the interferences are 
again negative above the A peak, and small (positive or 
negative) below the A peak. The overall increase of the 
cross section in comparison with the Delta pole contribu- 
tion is, as expected, much more significant than for the 
proton target. 

A feature of the HNV model is that it introduces the 
background not only for the pn^ and n7r+ final states, 
but also for the pn^ , that is for the isospin-3/2 channel. 
We observe, that in this channel the contribution of the 
background is at the level of 10%, which agrees with the 



result of [lj|. This justifies the neglect of the background 
in the isospin-3/2 channel as assumed in earlier works |7|, 
32| and explains why they were still successful 
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in describing the data. 



V. GENERAL FEATURES OF EFFECTIVE 
BACKGROUND 

As the next step, in Fig. [5] we present our results for the 
integrated cross section, with the kinematical cut for the 
nucleon-pion invariant mass W{Ntt) < 1.4 GeV, versus 
the neutrino energy. The calculations are made for vari- 
ous final states and compared with data and with some 



> 

O 

CM 

E 




1.1 1.2 1.3 1.' 
W, GeV 



1.5 



1.1 1.2 1.3 1.4 1.5 
W, GeV 



FIG. 4: Double differential cross section da /dEf_id cos 6^ for 
various final states for E^ = 1 GeV, cosO — 0.6. Contribu- 
tions of each diagram (a, c, e), as well as their coherent and 
incoherent sums (b, d, f) are shown. 
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previous theoretical results. 

The simplest channel to compare with is the scatter- 
ing on a proton target, because only one final state, p7r+, 
is possible. The full model calculation (solid line) ap- 
pears to be slightly above the Delta pole contribution 
(dash-dotted line) and coincides with the previous cal- 
culation [93 (dashed line labeled "Leitner 09") at small 
neutrino energies.^ With increasing E^^ the "Leitner 09" 
curve, as expected, increases more steeply than the full 
model curve because the calculation [9| was done with- 
out any kinematical cut, while our calculation implies 
WiN-K) < 1.4 GeV. Comparison with Fig. 5 in [14] 
shows, as expected, that the integrated cross sections is 
also very close to the original HNV result. 



FIG. 5: The integrated one pion cross section, with kinemati- 
cal cut VI^(A''7r) < 1.4 GeV versus neutrino energy for various 
final states. The full model calculations (solid line) are com- 
pared with the Delta pole contribution (dash-dotted line) and 
previous calculations of Leitner et al i9|] (dashed line) . 



For the reactions on the neutron two final states, nn'^ 
and p7r°, are possible. For both of them, the full model 
cross sections are close to the previous GiBUU results [9| , 
but have slightly different shapes. Keep in mind, that the 
W cuts are different. ^ 



^ This implementation is available in the current open-source ver- 
sion of GiBUU [ll. 



^ For the pir^ channel our result is also in agreement with the 
original HNV calculation [14l | (see Fig. 5 there), while for the 
n7r+ channel it is noticeably higher (at Ei, = 1.6 GeV our re- 
suh 0.12 ■ 10-38 cm2 versus HNV 0.08 • lO'^S cm^). To un- 
derstand this difference we compared our calculations for each 
diagram with the corresponding unpublished results of the HNV 
group. We found a very good agreement for all diagrams except 
cDp, which in our calculations appears to be around 1.7 times 



The ANL data can be described quite well for all 
channels and over the full energy range. This agree- 
ment is trivial for the pn^ channel because, as we men- 
tioned before, the axial form factors were fitted to them. 
The agreement for the two other channel shows that the 
model gives a very reasonable estimate for the integrated 
background. 

A few data points from the Gargamelle propane ex- 
periment at CERN PS 33| are also available for antineu- 
trino reactions on neutron Dn and nucleon PN, the latter 
being the sum over proton and neutron targets. Our 
full model calculations for the three possible final states, 
Dn — >■ fj,~^mr~, Dp -^ fi^wK^, and vp — >■ fi~^pTr~, are 
shown in Fig. |6] and compared to the Gargamelle exper- 
iment. The data on the pn^ channel are obtained as 

a{Dp — > /i^pvr^) = a{vN — )• fi^ N7r^)—a{Dn — )• fi^mr^) . 

The results of our calculations are very close to those 
presented in the HNV paper lj|. The cross section for 
mT~ is overestimated, while that for p7r~ shows a good 
agreement with the data. 

Comparing the curves for the full model and Delta pole 
contributions, one can define the effective background by 
subtraction: 

„eff-bgr „full model „Delta pole 

By definition, cr^n-bgr jj^dudes interference terms and 
thus can take positive or negative values. This effective 
background can be compared with the phenomenological 
background used in [9|. 

Fig. [7] shows the effective background for the three final 
states. For the p7r+ final state (long-dashed line), it is 
at the level of 0.04 • lO"'^^ cm^ for all values of E^. For 
a neutron target the background is growing with energy 
taking on about the same values for the pir'^ (solid line) 
and n7r+ (short-dashed line) final states. 

In earlier phenomenological approaches [7|, |9| , the as- 
sumption of the isospin-1/2 background was used, which 

presupposes iyJ^Z._i/2 = '^'^bgr-1/2'' ^^^^ curve is shown 
as a dash-dotted line. Thus, one can conclude, that the 
effective background does not follow the isospin-1/2 ap- 
proximation. Indeed, from the six diagrams (cDp, Np, 
cNp, CT, pp, pF) directly contributing to the back- 
ground, only one (Np) satisfies the isospin-1/2 hypoth- 
esis, that is, its Clebsch- Gordon coefficients are related 



smaller. Comparing the Clebsch— Gordon coefficients for various 
final states (see Table U} among themselves, one can notice that 
the cDp diagram contributes mainly to the n7r+ channel. Taking 
into account possible interferences, we attribute the difference in 
this channel to this contribution. 
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FIG. 6: The integrated one pion cross section, with kinemati- 
cal cut W{Ntt) < 1.4 GeV versus the antineutrino energy for 
various final states. 



as C ''+ = V2C ^. For others the corresponding rela- 
tion is different from \/2, and interference also plays an 
important role. 

Fig. [5] shows the ratio of the background to the full 
model cross section. To investigate how sensitive this 
result is to the W{N7t) cut, we plot two curves for each 
final state: with the cut WiNn) < 1.4 GeV (thin lines) 
and W{N'k) < 1.3 GeV (thick lines). As one can easily 
see, the results for these two cases are very close. For the 
pTT"*" channel, the ratio is large for low energy, but steeply 
falls down and does not exceed 10% for E^ > 1 GeV. For 
the neutron target the background is large and is at the 
level of 35% for the pn'^ channel and 50% for n7r+ one. 

Recall, that the background seen in neutrino reactions 
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FIG. 7: The effective background versus the neutrino energy 
for various final states. 
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FIG. 8: The ratios of the effective background to the full 
model cross section versus the neutrino energy for various 
final states for the cuts W{Ntt) < 1.4 GeV (thin lines) and 
W{N-!v) < 1.3 GeV (thick lines). 



(weak background) include vector, axial and vector-axial- 
interference parts. In electron reactions, only the vector 
part is present (electromagnetic background) and well 
constrained by data. In \^ it was assumed, that da^ 
and da^J^ have the same functional form as the vector 



part: 



dcTbgr 



d^bgr 



da. 



VA 



(1 



hN-K\^„V 



Ida,, 



(18) 



The vector parts were extracted from data independently 
for various channels. The coefficient b was the adjustable 
parameter fitted to the ANL data under the assumption 

Fig. [S] compares the effective background of the HNV 
model (solid lines) with the phenomenological back- 
ground P^ (dash-dotted lines) for the three final states. 



The agreement between the curves is reasonable for all fi- 
nal states up to E^, = 1.2 GeV, which justifies the model 
used in j9|]. 
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FIG. 9: The effective background versus the neutrino energy: 
comparison with the phenomenological calculations of [9|. 



VI. COMPARISON WITH THE ANL 
EXPERIMENT 

As we already mentioned, the data on neutrino-nucleon 
interactions are scarce and come from the late 1970s and 
early 1980s. In all these experiments wide band neutrino 
beams were incident on hydrogen and/or deuterium tar- 
gets. The most detailed sets of data are provided by the 
ANL 12-ft and the BNL 7-ft bubble chambers. The in- 
tegrated cross sections from ANL and BNL experiments 
were already used in the previous section. However, valu- 
able information comes also from the differential cross 



11 



sections. 

The deuteron effects for relevant neutrino energies were 
studied in J5|, wliere tliree different wave functions, cor- 
responding to the Hulthen, Bonn, and Paris NN models, 
were considered. It was shown, that these effects depend 
on the model used and generally suppress the cross sec- 
tion at Q2 < 0.1 GeV^ by no more than 8% for ANL 
experiment. At higher Q^ they are practically negligi- 
ble. For recent calculations including deuteron effects. 



see 



11 



34[ 



The GiBUU model implements the Argonne V18 NN 
potential. Within this model the deuteron effects are 
shown [20] to introduce only a minor correction to the 
ANL Q^ distribution — the correction is even smaller for 
BNL because of the higher neutrino energy. For invariant 
mass distributions, which are integrated over Q^, it would 
be even smaller. Thus, for the present calculations we 
neglect the effects of the deuteron structure. 



A. Transformation from events to absolute cross 
section 

Many reaction rates were presented in [2| not as ab- 
solute cross sections, but as events per some interval of 
the measured variable (Q^ or W ^ for example). This is 
mainly explained by the fact, that neutrino fluxes are 
not precisely known, but can only be determined with 
some accuracy which is hard to estimate. For the ANL 
experiment, for example, the flux is calculated from the 
measured multipion production cross section on a beryl- 
lium target and is given in [35[. It is clear, however, 
that the transformation coefficient fc from the number 
of events per unit energy to the absolute cross section 
is unique for a given experiment and (for a perfect ex- 
periment) must be the same for each reaction channel 
ch — {pn'^ ,pit'^ ,mT~^) for all distributions. 

The ANL experiment provides data for the distribution 
of observed (also called raw) events Nch{Eu) in neutrino 
energy |^ (see Fig. 7 there), "^ which is 

fch ■ N,h{E,) = (T,„t(,;,) {E,) ■ flux{E,) ■ k . (19) 

The rate correction coefficients fch account for experi- 
mental backgrounds and losses. They can be extracted 
from the summary of rate corrections [^ (see Table I 



^ Note, that neither ANL nor BNL provide distributions corrected 
for the experimental backgrounds, which, ideally, would be the 
subject of comparison with the theory. 



there) or calculated as ratios of corrected to raw events 
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JvTT« - 2023 ~ -L-04y, /„^+ — T^TTTT^ - i./4i 



206.2 



(20) 

With this data in hand and considering the flux and cross 
section as experimentally determined, we are able to cal- 
culate the coefficient k for different final states, compare 
the results and thus estimate the accuracy intrinsic to 
the experiment and its consistency. 

Fig. [ini shows the coefficients k determined from 
Eq. (fT9|) and the ANL data. The flux is taken from the 
histogram in [35J. The experimental points for the cross 
sections are interpolated with splines. The errors for the 
cross sections are used to estimate the errorbands of the 
coefficient. 

For each channel the middle curve shows the central 
value of k which corresponds to the central value of the 
cross section. The lower (upper) curves correspond to 
the maximal (minimal) values of the cross section and 
serve as error bands for k. For each channel k is fltted as 
a constant value, each point weighted with its maximal 
error. The results 

fcP"" -(370 ±311 ^"^^^^/"-^ 
^r.^.-(380±31)^S^^ 



(21) 



fcr.^,^ (380 ±31)^^5^ 

10 ■^'^ CTO^ 



are shown as straight lines in Fig. [101 The factor 0.1 
comes from the E^, binning in Fig. 1101 

Comparing kANL with the corresponding values for 
each channel separately, we conclude that an accuracy 
of around (453 - 370) * 2/(453 + 370) « 20% should be 
attributed to it. This value is consistent with the ANL 
flux uncertainty estimated in [2] as 15%. 

The procedure considered here is inverse to what ex- 
perimentalists do to determine the absolute cross sec- 
tions. Here its main purpose is to estimate a reasonable 
accuracy requirement for fitting theoretical curves to the 
data. We conclude, that for the ANL experiment an 
agreement within 20% should be considered as perfect. 



B. Q^ distribution 

For the pir^ channel the Q^-distribution is given by 
the ANL experiment as an absolute cross section dcr/dQ"^ 
for events with the invariant nucleon-pion mass cut 
W^ < 1.4 GeV and neutrino energy cut 0.5 GeV < E^, < 
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FIG. 10: The ANL transformation coefficient, determined 
from the data on neutrino event distributions for various final 
states. The three lines reflect an error band obtained from 
the errors of the published data. 
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FIG. 11: Cross section da/dQ^ averaged over the ANL neu- 
trino energy flux for the final state ^Tpir'^ . The integration 
over W is performed with a VK < 1.4 GeV cut, in agreement 
with the experiment [2|. Data are shown as filled squares. 



the experimental histograms. The Delta pole contribu- 
tions to the pTT^ and m:'^ channels are noticeably below 
the data. The extra contribution from the background 
adds around 50% for the pn'^ channel, which overshoots 
the data at low Q^ . The general agreement of our curve 
with the data is very good. For the n7r+ channel, the 
background contribution adds 100% to the cross section, 
which, however, is still not enough to reach the experi- 
mentally observed values at low Q^ . This could hint at 
a contribution of the higher mass isospin-1/2 resonances, 
not considered in this work, which may decay into one 
pion final state. With the estimated intrinsic uncertainty 
of 20%, the overall agreement should be considered as 
good. 



W distributions 



6 GeV. Our results, presented in Fig. [TTl show a good 
agreement with the experimental data. 

The same distribution, but without energy cut, as well 
as those for the pi:^ and ni:^ channels, are presented as 
events per Q^ interval. As mentioned before, the axial 
form factors of the theoretical model were fitted to the 



pTT 



channel of the ANL data; we thus use kANL = k^i\ 



ANL 



given in Eq. (PT|) as our transformation coefficient for all 
distributions. In comparing our theoretical results with 
the data we therefore normalize them to each other by 
multiplying the theory results with the factor kANLJ fch 
with fch for the various channels given in Eq. pO]). 
Our results are shown in Fig. [T^] and compared with 



Next, we present data for invariant mass distributions. 
In previous theoretical investigations [4|, |6|, |lJ] , only dis- 
tributions versus pion-nucleon invariant mass W{N'k) 
were calculated. The ANL and BNL experimental data 
are available also for nucleon-muon W{plN) and pion- 
muon iy(/i7r) combinations. These additional data can 
be used to constrain the theory even further. 

Fig. [T3k shows the W{Ntt) distribution for the p7r+ 
channel. The agreement of our full model calculations 
(solid curve) with the histogram is very good. Of inter- 
est is the region of low W , near the one pion production 
threshold. In this region, the experimental data show a 
noticeable rise with increasing W ^ which is in agreement 
with the full model prediction. The Delta pole contribu- 
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FIG. 12: Cross section da /dQ^ averaged over the ANL neu- 
trino energy flux for the final states: (a) fj,~p-K'^ , (b) ix'pn'^, 
and (c) n~mT^ . The integration over W is performed with 
a W < 1.4 GeV cut, in agreement with the experiment [2|. 
Data are shown as histograms. 



FIG. 13: The nucleon-pion invariant mass distributions, av- 
eraged over the ANL flux. The fuU model calculations (solid 
curve) and Delta pole contribution (dash-dotted curve) are 
shown. The experimental data from [2] are shown as his- 
tograms. 



tion (dash-dotted curve) , on the other hand, grows rather 
slovifly. 

For the p7r° and nn^ final states, as shown in 
Figs. [13b.c. the agreement of the full model with the 
histogram is reasonable. For the pir'^ channel the full 
model overestimates events in the Delta peak region and 
underestimates them immediately above this peak. For 
the nTT"*" channel, the data are underestimated below the 
Delta peak. 

At low W the background gives a noticeable contribu- 



tion for both channels, in line with the data. For differ- 
ent final states the background contributions above the 
Delta peak are very different: small negative for p7r+, 
very small for pTr" and positive for nir^ . 

While the M^(iV7r) distributions are mainly sensitive to 
the A excitation, the distributions W{fiN) and W{^tt) 
test the angular distribution of the vN interaction. 

The W{iJ,N) and W{fnr) distributions shown in 
Figs. [131 [m also agree reasonably well with our calcu- 
lations. Recall, that in the ANL experiment the pTr"*" 
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FIG. 14: The muon-nucleon invariant mass distributions, av- 
eraged over the ANL flux. The full model calculations (solid 
curve) and Delta pole contribution (dash-dotted curve) are 
shown. The experimental data from 2] are shown as his- 
tograms. 



FIG. 15: The muon-pion invariant mass distributions, aver- 
aged over the ANL flux. The full model calculations (solid 
curve) and Delta pole contribution (dash-dotted curve) are 
shown. The experimental data from [2] are shown as his- 
tograms. 



data are presented for the whole neutrino energy flux, 
which only vanishes at E^ = 6 GeV. This can explain 
the large tail in this distribution. For the pn^ and mr^ 
final states, on the other hand, the experimental data (as 
well as our calculations) are limited to E^ < 1.5 GeV, so 
that the large W{iiN) are not kinematically accessible. 

The full model and the Delta pole terms give curves 
of similar form, but different magnitude in the various 
isospin channels. For the pn'^ channel the agreement of 
our the full model with the data is very good. In the pyr" 



channel the data are underestimated at low W{fiN). For 
the mr~^ channel the data are underestimated at both low 
WifiN) and low Wifin). 

The overall results clearly indicate the necessity to in- 
clude the nonresonant contribution is addition to that 
of the Delta pole. Even with some underestimation, the 
full model curves show much better agreement with the 
whole set of data than the Delta pole terms alone. The 
background adds around 10% to the Dp cross section for 
the p7r+ channel, around 50% for the pir'^ one and around 



15 



100% for the mT+ one. 

From Figs. [T5HT51 we conclude, that within the exper- 
imental accuracy available, the data presented can dis- 
criminate between the Delta pole and the full model 
curves and are compatible with the HNV background 
model. 



VII. COMPARISON WITH THE BNL 
EXPERIMENT 

A. Transformation from events to absolute cross 
section 

The results of the BNL experiment on Q^ or W distri- 
butions, like those of the ANL experiment, are presented 
as events. Data on the distribution of events in neutrino 
energy are also available [3[ (see Fig. 2 there), so we use 
the same procedure as before to determine the transfor- 
mation coefficient and estimate its accuracy. 

The rate correction coefficients fch for BNL are given 
in [il (see Table II there) 



/, 



BNL 



1.12, L^J^ = 1.05, /, 



BNL 

+ 



0.89 . (22) 



Fig. [TB] shows the coefficient k determined from 
Eq. (IT91) and the BNL data. The flux is calculated from 
the observed quasielastic events and given in 36|. The 
experimental points for the cross sections are interpo- 
lated with splines. The error bars of the cross sections 
are used to estimate the error bands of the coefficient. 

For each channel, the middle curve with data points on 
it shows the central value of fc, while the lower (upper) 
curves serve as error bands for k. For each channel, k is 
fitted as a constant value, each point weighted with its 
maximal error. The results 
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FIG. 16: The BNL transformation coefficient, determined 
from the data on neutrino event distributions for various final 
states. The three lines reflect an error band obtained from 
the errors of the published data. 
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(23) 



are shown as straight lines in Fig. [161 

Comparing ksNL for the three channels, we attribute 
an accuracy of around (544 - 518) * 2/(544 + 518) = 
5% to ksNL- We conclude, that when comparing with 
the BNL experiment, an agreement within 5% should be 
considered as perfect. 



B. Q^ distribution 



The data on the Q^ distribution for the pir^ chan- 
nel are presented in [3| for the neutrino energy cut 
0.5 GeV < i?^ < 6 GeV and the invariant mass cut 
W < 1.4 GeV. Since the latter corresponds to the range 
of applicability of the HNV model, we can normalize the 
area under the full model theoretical curve to that under 
the experimental data. This is another way to estimate 
the transformation coefficient fc^^vi [which for a perfect 
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experiment must be equal to those in Eq. (l23l) ]: 
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(24) 



This value exceeds the typical value of fcsATL determined 
in Eq. ((23)) by more than 30% (even more for the p7r+ 
channel), and thus cannot be considered as consistent. 

Speculating about the possible origin of this inconsis- 
tency, we note that the values (|23|) for different channels 
agree among themselves quite well, which may hint at 
the consistent treatment of the rate correction coefficients 
fch- Keeping also in mind, that the flux is the same for 
various channels, a possible way to explain the above in- 
consistency would be to suppose that the cross sections 
are overestimated by 30%. This is exactly the difference 
between ANL and BNL integrated cross sections. Thus, 
by reducing crf3^i^-. by 30% one would simultaneously 
reach agreement with ANL and obtain k'j^j^j^ consistent 
with k^^. 

With the data as they are, we reestimate a realistic 
uncertainty as 30%. 

Thus, for the BNL experiment we aim at shape-only 
comparison and hereinafter use fc^JVL p4|) as our trans- 
formation coefficient. 

The comparison of the Q^ distribution with our cal- 
culations is shown in Fig. [T71 For the p7r+ channel, 
as expected, the background gives only a small contri- 
bution in addition to the leading Delta pole term. In 
agreement with earlier calculations |32j , the peak of the 
curve is located in the region Q^ = 0.1 — 0.12 GeV^ 
and is shifted with respect to the data peak at Q^ = 
0.18 — 0.2 GeV^. This disagreement has been known for 
a long time with similar results obtained within various 
models [3|, l20|, 132, |37| and is not resolved. Inclusion of 
the background does not change the peak position. 

The data for the other channels are given in [3| without 
any cut on W, and for the whole neutrino flux (that is 
0.34 GeV < E^ <6 GeV) as events per 0.1 GeV^ inter- 
val. They are compared with our calculations in Fig. I18[ 
where the same transformation coefficient k^^j^/f^^f^^ 
[Eq. p4)) ] is used. The same coefficient is also used fur- 
ther for various W distributions. 

For the p7r+ channel, as expected, our full model calcu- 
lations (solid line) is below the experimental histogram. 
The area under our curve is 15% below the area un- 
der the histogram. Since the data include higher invari- 
ant masses whereas our calculations contain a cutoff of 
1.4 GeV, corresponding to the range of validity of the 
HNV model, this implies that 15% of all events in this 
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FIG. 17: Cross section da/dQ^ averaged over the BNL neu- 
trino energy flux for the final state fi'pn'^. The integration 
is performed with the W{Nn) < 1.4 GeV cut, in agreement 
with the experimental data [3], which are shown as histogram. 



channel should be attributed to the higher mass isospin- 
3/2 resonances, such as P33(1600),S'3i(1620),D33(1700), 
and their interferences with the background, which are 
not considered here. 

For the pir'^ and mr^ channel the background signifi- 
cantly increases the cross section in comparison with the 
Delta pole contribution. However, our full model curves 
are still much lower than the histograms, which indicates 
a large contribution of higher mass isospin-1/2 and -3/2 
resonances and their interferences. The relative impor- 
tance of these events is estimated by comparing the areas 
under the theoretical curve and experimental histogram, 
as it was described for the pTr"*" channel, and it appears 
to be 43% for pn^ and 46% for mr~^. This will also be 
demonstrated further in the WiNir) invariant mass dis- 
tribution. 

Notice also, that in the p7r° channel [Fig. [TST b)] the 
peak of our curve is shifted to the left with respect to 
the histogram. This effect is the same as in pn^ channel 
in Fig. [T71 but it is revealed here with less significance 
because of the larger Q^ binning. 

C. W distribution 

Now we proceed with calculating the invariant mass 
distributions. Fig. [T9l shows the W{NTr) distribution for 
the three final states. Our calculations are done only up 
to Vt^ < 1.4 GeV, which is the range of applicability of the 
HNV model, while the experimental data are available 
also for higher W. 

For the pTr"*" channel, our calculations show a very good 
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FIG. 18: The da jdQ^ cross section averaged over the BNL 
neutrino energy flux for the final states: (a) jj,~ptt'^, (b) 
^Tp-K^, and (c) fi~mv'^. The integration is performed with 
the W{N-k) < 1.4 GeV cut, corresponding to the range of 
apphcabihty of the HNV model . The experimental data [Sj] 
shown as histograms are without W cut. 
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FIG. 19: The nucleon-pion invariant mass distributions, av- 
eraged over the BNL flux. The full model calculations (solid 
curve) and Delta pole contribution (dash-dotted curve) are 
shown with the cut W{NTr) < 1.4 GeV. The experimental 
BNL data [3] which do not contain this cut are shown as his- 
tograms. 



agreement with the data. In the region 1.05 GeV < W < 
1.4 GeV the area under our curve coincides with the area 
under the histogram with an accuracy better than 1%. 

To estimate how many events belong to the high 
W{NTr) region, we calculate the area under the histogram 
with the cut W{Ntt) < 2.0 GeV and compare it to that 
with the cut W{Ntt) < 1.4 GeV. We find that > 10% of 
all events belong to the H^(iV7r) > 1.4 GeV region. This 
is in agreement with the conclusion previously derived 



from the Q^ distributions. 

For both the p7r° and mr^ channels, the background is 
essential, especially at very low W{NTr), where it signifi- 
cantly increases the cross sections. The overall agreement 
of our calculations with the data is reasonable. Similar 
to our comparison with the ANL experiment, for the pn'^ 
channel the full model overestimates the data in the re- 
gion of Delta peak and underestimates them immediately 
above this peak. For the nTr"*" channel, the data are un- 
derestimated below the Delta peak. The percentage of 
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events with T^(iV7r) > 1.4 GeV is estimated to be 45% 
for the pTr" and 44% for the n7r+, which is in good agree- 
ment with the previous results obtained when discussing 
the Q^ distribution. 

As in our calculations for the ANL experiment, in the 
W region above the Delta peak the background contribu- 
tion is different for different channels: negative for p7r+, 
small for yn^ and positive for n7r+. The comparison of 
the two latter channels shows that the effective back- 
ground da^^"^ jdW^Ni^) as it is described within the HNV 
model does not support the isospin-1/2 hypothesis. 



Delta pole 

full model 

BNL 




1.5 2 

W(nN), GeV 



FIG. 20: The muon-nucleon invariant mass distributions, av- 
eraged over the BNL flux, for the events with WiNix") < 
1.4 GeV compared to the BNL data fSl] shown as histograms. 

Fig. [50] shows the muon-nucleon invariant mass distri- 



bution compared with the BNL data subject to selection 
WiNir) < 1.4 GeV. 

For the pn^ channel, the background slightly increases 
the cross section. For the pn^ and mr^ channels it 
gives a significant contribution, which noticeably im- 
proves agreement with the histograms. For the pir^ chan- 
nel, the data show an excess of events over our curve at 
low W{iiN). As explained in [3], "this excess comes in 
part from the misidcntified lyn — >■ v'K~p and nn — >■ np-K~ 
events which belong to the experimental background" not 
subtracted from the data. 

This generally means that the experimental back- 
ground is nonuniform and requires a more detailed treat- 
ment. Indeed, all experimental histograms show the dis- 
tributions of observed (raw) events. The rate correction 
factors f^h^^ are provided experimentally as constant 
factors that do not depend on kinematics. In this way, we 
correct the total number of events, but the nonuniformity 
of the experimental corrections remains unaccounted for 
and can reveal itself in all distributions. 

For the uTr"*" channel, the agreement is good in the 
region of high W{iJ,N), but again we underestimate the 
data for low W{^iN). 

Fig. [21] shows the muon-pion invariant mass distribu- 
tion compared with the BNL data subject to selection 
W{Ntt) < 1.4 GeV. The agreement of our calculations 
with the data is good for the pir^ and pir^ channels, while 
for the nTT"*" one we underestimate the data at low W{fiTT). 

All the three W distributions taken together clearly 
show the importance of background terms in reaching a 
reasonable description of the data for all the three in- 
variant mass distributions considered. Delta pole term 
alone, showing a good agreement with the data for the 
p7r+ channel, significantly underestimates all the data for 
the pir^ and nTr"*" channels. 

Some discrepancies between the full model calculations 
and data at low W^s [for aU three W^(iV7r), W{iiN) and 
W{fiTT)] for both ANL and BNL may point to the neces- 
sity to improve the model. As discussed earlier in this 
section, the BNL experimentalists themselves attributed 
the excess of events in W{fiN) distribution of the j>7r° 
channel to some misidentified events. Similar excess in 
other channels for low W may have a similar origin. 



VIII. DISCUSSION AND CONCLUSION 

In Sec. [nj we have outlined the HNV model and the 
phenomenological form factors for various diagrams. As 
ensured by comparison with electroproduction in Sec. IIIII 
the model is applicable up to nucleon-pion invariant mass 
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FIG. 21: The muon-pion invariant mass distributions, av- 
eraged over the BNL flux, for the events with W{Nn) < 
1.4 GeV compared to the BNL data [3i] shown as histograms. 



of W{Nn) < 1.4 GeV and provides about the same level 
of accuracy as the MAID model. 

For neutrino reactions, as discussed in Sec. IIVI the 
model predicts a background contribution, which is small 
for the p7r+ channel (at the level below 10% for neutrino 
energies above 1 GeV), at the level of 30% (with respect 
to the full model cross section, that is 50% with respect 
to the Delta pole contribution) for the pTr" channel and 
at the level of 50% (with respect to the full model cross 
section, that is 100% with respect to the Delta pole con- 
tribution) for the nTT"*" channel. The effective background 
does not satisfy the isospin-1/2 hypothesis. 



The HNV model describes the available data set on 
neutrino and antineutrino reactions on nucleons reason- 
ably well, with an accuracy that approximately corre- 
sponds to the accuracy with which different data agree 
among themselves. 

For neutrinos the absolute values of the integrated 
cross section are available from the ANL and BNL exper- 
iment, with the BNL data being systematically higher. 



The Delta axial form factors were fitted in lj| to the 
ANL integrated cross section and the Q^ distribution for 
the p7r+ channel so that the full model agrees with them 
by definition. Agreement with the ANL data for the pn'-' 
and m:'^ channels shows that the model gives a good de- 
scription of the background. The BNL data lie around 
30% higher for all channels. 

For antineutrinos only very few data are available from 
the Gargamelle experiments. The agreement of the full 
model calculations with the data is good for the pTT~ 
channel and overestimates the data for the rnr^ one. 

For most of the differential cross sections available ex- 
perimentally no information about the absolute value of 
the cross section is available; the data are presented as 
raw events per Q^ or W interval. By estimating the 
transformation coefficients for the ANL and BNL exper- 
iments in Sees. iVIblVIIb . we have shown that the real- 
istic accuracy for the ANL experiment is around 20%, 
with the uncertainties coming from the inconsistencies in 
the various channels. For the BNL experiment the vari- 
ous channels agree within 5%; here the overall accuracy 
about 30% comes from the disagreement of the Q^ distri- 
bution, discussed in Sec. iVIIb . with the integrated event 
distribution. 



As discussed in Sees. IVIIVIIl the overall agreement 
with the data is perfect for the pn'^ channel and reason- 
able for the p7r° and nn^ ones. For the latter two chan- 
nels the full model calculations systematically underes- 
timate the experimental histograms at low W^s. These 
discrepancies may hint at the necessity to improve the 
model, but they may as well come from the recognized 
nonuniformity of the experimental background. 

Even without absolute normalization, and even tak- 
ing into account the discrepancies observed, all the ANL 
and BNL data taken together are able to discriminate 
between the full model and the leading Delta pole con- 
tribution. When all three final states accessible for neu- 
trino reactions {pir^ , pir^, mr^) are considered, the data 
definitely demand the nonresonant background and favor 
the full model calculations. 
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